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ABSTRACT 

Experinental transient electromagnetic fields measurements have been 
performed on a Transall C160 aircraft during in-flight lightning strikes. The 
data allow a test of the predictive capabilities of a three dimensional 
time-domain finite-difference code (ALICE) developed at ONERA in order to 
investigate lightning aircraft interactions. Using a transfer function 
technique in the 3D code it is shown that a bi-leader attached to an aircraft 
can be simulated by a linear model and so electromagnetic fields can be 
calculated everywhere on the vehicle. Comparison of experimental and 
numerical results have been made for several lightning strikes. Skin current 
density and electromagnetic fields distributions are discussed in detail. 


1. INTRODUCTION 

The electromagnetic pulses associated with a lightning threat can induce 
large charges and currents on the external skin of an aircraft in-flight 
that, in turn, couple into the interior of the vehicle where they can cause 
failures or temporary upsets in mission critical subsystems, if these 
subsystems are not sufficiently hardened. The recent introduction of advanced 
composite materials, sophisticated low level flights and mission critical 
electronics into aerospace vehicles has raised concerns that these may 
increase the vehicle's inherent susceptibility/vulnerability to transient 
electromagnetic fields. 

Because of this concern, DRET (Direction des Recherches, Etudes et 
Techniques) has been sponsoring, for several years, a research program to 
investigate aircraft interaction with direct lightning strikes and to develop 
analytical and numerical tools necessary for a theoretical determination or 
prediction of lightning electromagnetic interactions with a vehicle. 

The most recent in-flight program was conducted in France with an 
instrumented C160 Transall aircraft [1]. 

The data obtained allowed a test of the predictive capabilities of a 
three-dimensional finite-difference computer code set up at ONERA in order to 
calculate charge and current distributions on the external skin of an 
in-flight vehicle struck by lightning. 

Magnetic field coupling through a carbon composite panel located on the 
fuselage has also been studied [2]. The experimental and numerical results 
are compared for several strikes. 
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2. IN-FLIGHT INSTRUMENTATION DESCRIPTION 


During Summer 1988, ONERA has performed a complete airborne lightning 
characterization program, including sensors design, instrumentation lay out, 
extensive ground test of the whole aircraft set in a coaxial return path and 
in-flight measurements. The sensors were implemented on the aircraft as shown 
in figure 1. Four resistive current shunts (I) were mounted on the aircraft. 
One was fixed to the nose-boom, two were fixed to the wing tips and one was 
located at the base of the tail-boom. Two derivative current sensors (I-Dot) 
were also installed in the middle of the front and tail-booms. Seven 
locations were instrumented with electromagnetic sensors for E, H, E-Dot, 
H-Dot external fields measurements (noted EH1 to EH7 in figure 1). An 
internal magnetic field measurement has been achieved at the center of the 
fuselage, behind a carbon composite panel located on the left side of the 
aircraft. Seven video cameras were used to record the scene all around the 
aircraft. One of then was a 200 frames/second camera located under the right 
wing. A more detailed description of the instrumentation and of the sensors 
characteristics is given in [1]. 

3. NUMERICAL SIMULATION METHOD 

The numerical code "ALICE" employed in this study can take into account 
lightning electromagnetic interactions with a vehicle (nearby and direct 
strikes, electromagnetic coupling through lossy structures, coupling through 
resistive and isotropic joints, ...). 

3.1. FINITE-DIFFERENCE TIME-DOMAIN METHOD 


The code is based on an explicite Finite-Difference Time-Domain (FD-TD) 
algorithm developed by K.S. Yee [3] in cartesian coordinates. The goal of the 
FD-TD method is to model the propagation of an electromagnetic wave into a 
volume of space containing a dielectric or conducting structure. The 
procedure produces a finite difference approximation to the solution of 
Maxwell's differential curl equations implemented on rectangular cross 
section-unit-cell space lattice. 

The differential time domain form of Maxwell's equations in an isotropic 
medium is: 

dB 


V x E = 

(1) 

at 


dD 

(2) 

7 x H = J + — 


dt 

with: B = i* H ; D = l E. 

O O 


These equations are written in the code in a finite difference form. For a 
parallelepipedic .cell space lattice, this procedure involves positioning the 
components of E and H about a unit cell of the lattice, as shown in figure 2, 
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and evaluating E and H at alternate half time steps. In this manner, centered 
difference expressions can be used for both the space and time derivatives to 
attain second-order accuracy in the space and time increments. The size of 
the unit cell determines the time step At which, for numerical stability 
condition, must satisfy the criterion: 


At 2 < 


I ( » 

o o 


1/Ax 2 + 1/Ay 2 + 1/Az 2 


(3) 


where Ax, Ay and Az are the dimensions of a cell (figure 2). 


3.2. LATTICE TRUNCATION CONDITIONS 

To minimize computer storage, the lattice must be truncated as close as 
possible to the modelled structure. However, care must be exercised because 
this condition must not cause excessive spurious reflection of waves 
scattered outward by the structure. 

The field components at the lattice truncation planes must be computed using 
an auxiliary radiative truncation condition. The method used in "ALICE" is 
achieved by using highly absorbing boundary conditions to simulate free 
space. This approximation is successful when the lattice is at least twice 
the largest vehicle dimension in each direction of space. 

3.3. THIN WIRE FORMALISM 


FD-TD codes cannot resolve coupling to wires of diameters less than a mesh 
size unless special features are added. The thin-strut formalism developed by 
R. Holland [4] allows such resolution and has been introduced in "ALICE". On 
a wire segment parallel to the x-axis (X, ^ coordinate location of the 
wire - figure 2) . The current I and the charge Q per unit length are 
calculated by finite differencing the two following equations: 


dQ ai 

3t dx 

ai , as <V 

at dx l 

> 

where 


(4) 

(5) 


C is the light velocity 

<E x > is the electric field component averaged over the cross section of the 
cell surrounding the wire. 

L^ is the in-cell inductance, the inductance per unit length a thin wire 
would have with respect to an enclosing conductor half a mesh removed. 


Thus, a cell with a wire running through it requires eight quantities (six 
fields components, the current 1 and the charge density Q) to be advanced 
each program cycle. 
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3.4. TRANSALL AIRCRAFT MODELLING 


The Transall aircraft was modelled using three-dimensional rectangular cells 
having dimensions of Ax = 45 cm, Ay = 90 cm and Az = 45 cm that were part of 
a total cell space of 121 x 70 x 79 cells. The relatively small cell size 
allows a reasonable condition of the aircraft geometry to be made as figure 3 
reveals. The time step At employed is found from the condition (3) and has 
been chosen equal to 0.9 ns. For the external fields calculations a uniform 
mesh has been used inside a box that just encloses the modelled 

aircraft-extremities and slowly varying around it (adjacent cells differed in 
size by 20%). For the internal fields calculations, a varying mesh has been 
chosen on the aircraft in order to properly take into account the right 
composite door dimensions and a second one in space around the aircraft as in 
the previous meshing. 

The metallic parts of the vehicle are modelled by the necessary condition: 

n A E = 0 (6) 


on the structure, n is a unit vector perpendicular to the surface. 

The connecting lightning channels are modelled using thin wire formalism 
describe in 5 3.3. 

A test of the predictive capabilities of the computer code "ALICE” has been 
performed, several years ago, by comparing experimental and numerical results 
obtained on an aircraft mock-up [5]. 

4. COMPARISON OF NUMERICAL PREDICTIONS AND IN-FLIGHT MEASUREMENTS 

We give here the comparison of experimental and theoretical results obtained 
for a nose-wing lightning strike obtained during the 1988 in-flight 
experiment. For this event electromagnetic signatures on the sensors and 
photographic evidence from the on-board video cameras indicated that the 
aircraft initiated a bi-leader mechanism. The sequence of the phenomena is 
the following (figure 4a) : 

- at the beginning (between A and B) the aircraft is immersed in an ambiant 
electrostatic field, 

- at the nose-boom, the electric field is high enough to initiate a positive 
leader which propagates towards ground (between B and C), 

- at the right wing tip, a negative leader is initiated and propagates, by 
steps, towards the rear part of the fuselage (between C and E) . The 
direction of propagation is parallel to the fuselage as shown in figure 3. 
It is to be noticed that the negative leader propagation generates current 
pulses which are superimposed on a slowly varying current component 
(figure 4b) . 

The numerical simulation of this event has been achieved for the phenomena 
occuring between points C and E in figure 4a. During this phase, the positive 
leader could be considered as having an infinite length and the negative 
leader could be considered as having a 30 m characteristic length 
corresponding to its first step [6]. 
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The numerical calculations have been undertaken in two steps using a 
frequency domain transfer function technique. 

4.1. ELEMENTARY PULSES MODELLING 

The principle of the transfer function method is the following: 

a) We assume the lightning aircraft interaction is a linear mechanism. So, 
for a given geometry of the vehicle attached to lightning channels (as 
shown in figure 3) one can put, in the 3-D code, a current generator in 
place of the negative leader tip and which injects a known current pulse 

I (t) towards the aircraft. In one point p of the vehicle (taken as a 
0 © 
reference point) a magnetic field response H o (t) can be calculated and 

using a FFT of H (t) and I (t) one can then determine (in the frequency 
0 0 

domain) a transfer function T(F) between these two points: 

T (F) = H (F)/I (F) (7) 

0 0 

where F is the frequency and H (F) , I (F) are the FFT of H (t), I (t) 

0 0 o o 

respectively. 

We generally use a standard sine-squared pulse for I 0 (t): 

I (t) = I sin 2 (nt/2T) (8) 

O M 

where t is chosen in function of the frequency content of the in-flight 
measured responses. 

b) If we consider the response H a ( t ) which has been obtained, at the same 
point P # on the aircraft, during in-flight experiment, its FFT H 3 (F) 
allows the calculation, in the frequency domain, of an equivalent current 
generator 1(F) located at the negative leader extremity and defined by: 

1(F) = H s (F)/T(F) (9) 

and in time domain: 

I (t) = I FT [1(F)] (10) 

c) This current waveform I(t), simulating the negative leader mechanism, is 
introduced in the 3-D code and then the temporal responses at every 
instrumented points on the aircraft (figure 1) can be computed and 
compared to experimental data. 

This technique is illustrated in figure 5 to 8, for the measurement 
current pulse noted "pulse n'7" in figure 4b. The pulse noted "pulse n*2" 
has been studied previously [5] using the same technique. 

Figure 5 represents the comparison between magnetic field pulses measured and 
calculated on the left side (curves a) , on the right side (curves b) and on 
the top (curves c) of the forward fuselage. 

Figure 6 represents the comparison between magnetic field pulses measured and 
calculated inside the fuselage, behind the carbon composite door. 
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Figure 7 represents the comparison between electric field pulses measured and 
calculated on the top of the forward fuselage. 

Figure 8 represents the comparison between the current pulses measured and 
calculated on the right wing tip (curves a) and on the nose boom (curves b) . 
The magnetic field H 3 measured on the top of the fuselage has been used as 
reference to compute the transfer function T(F). The comparisons show a good 
agreement between measured and calculated waveforms at different locations on 
the aircraft, the magnitude discrepancy is better than 3 dB for all the 
fields and currents. 

4.2. PULSE SERIES MODELLING 


Knowing the entire strike response of a sensor in one point on the aircraft, 
we can calculate the responses everywhere on the aircraft using the transfer 
function technique described previously. 

In the phase a) of S 4.1 we can calculate the transient response H (t) at the 
reference point P but also the transient responses H (t) at different points 
P of interest on the aircraft. 

We can then determine (in the frequency domain) the transfer functions T (F) 
beetween the points P and the reference point P , such as: 

i 0 

T, (F) = H <F)/H (F) (11) 

1 1 O 

where H t (F) and H o (F) are the FFT of H j (t) and ( t ) respectively. 

If we consider the pulses series H a (t) (figure 9) which has been measured 
during in-flight experiment at point P^ on the top of the front fuselage 
(figure 1), its FFT H (F) allows the calculation, in the frequency domain, of 
the responses H t (F) at each point : 

H ( (F) = Tj (F) .H 3 (F) (12) 

This technique is illustrated in figures 10 to 12. 

Figure 10 represents the comparison between the magnetic fields measured 
(curve a) and calculated (curve b) on the left side of the forward fuselage. 

Figure 11 represents the comparison between the magnetic fields measured 
(curve a) and calculated (curve b) inside the fuselage, behind the carbon 
composite door. The pulses of the measured field (curve a) are superimposed 
on a CW component due to an induced effect on the sensor of the aircraft 
power supply. 

Figure 12 represents the comparison between the currents measured (curve a) 
and calculated (curve b) at the right wing tip. As explained in S 4, the 
measured current pulses are superimposed on a slowly varying current 
component. This low frequency component has not been measured by the magnetic 
field sensor H (cut-off frequency of 100 Hz - figure 9) and so can not be 

3 
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calculated (figure 12b) using relation (12) . 

5. CONCLUSION 

The goal of this paper is to present a 3-D FD-TD computer code developed at 
ONERA to model ligthning interaction with an aircraft. The predictive 
capabilities of the code have been quantified, a few years ago, by comparing 
experimental and computed results obtained on an aircraft mock-up and on an 
helicopter carbon composite fuselage. 

We show in this paper that the code may be a powerful tool to analyse 
in-flight data. Using a transfer function technique, it has been shown that a 
bi-leader attached to an aircraft can be simulated by a linear model and so 
E.M. fields can be calculated everywhere on the vehicle. 

Furthermore, knowing the entire strike response of a sensor in one point on 
the aircraft, we can calculate the responses everywhere on the aircraft. 
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a) On the left side of the forward 
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b) On the right side of the forward 
fuselage. 


*1 H3 



c) On the top of the forward fuselage 

Fig. 5 : Comparison of measured and 
calculated magnetic field 
pulses: 



Fig. 6 : Comparison of the magnetic 
field pulses measured and 
calculated inside the fuse- 
lage behind the carbon com- 
posite door. 
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